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ABSTRACT 
About 900 fast radio burst (FRB) sources have been detected till now, among which 29 FRBs are 


found to burst out repeatedly. Although a firm connection between at least some FRBs and magnetars 


has been established, the trigger mechanism and radiation process in these enigma phenomena are 
still highly controversial. In this study, we build a sample of 16 repeating FRBs from which at 
least five bursts have been detected, including the most active four repeaters of FRBs 20121102A, 
20180916B, 20190520B and 20201124A. Various key parameters of their bursts are collected from the 
literature, which include the trigger time, pulse width, dispersion measure (DM), Faraday rotation 
measure (RM), bandwidth, waiting time, peak flux and fluence. The distribution and time evolution 
of these parameters are investigated. Potential correlations between various parameter pairs are also 
extensively explored. The behaviors of different repeaters are then compared. It is found that the 
DM of FRB 20121102A seems to increase continuously on a long timescale. While the DM of most 
repeaters varies in a narrow range of +3 cm ^? pc, FRB 20190520B is found to have a large variation 
range of +12 cm ? pc. A linear correlation between DM and the star formation rate is established. 
RM evolves with time in a much more chaos behavior in different repeaters. Generally, the waiting 
time shows a bimodal distribution in each source. The implications of these features to the underlying 


physics are discussed. 


Keywords: Radio transient sources (2008); Astrostatistics (1882); Model selection (1912); High energy 
astrophysics (739); Burst astrophysics (187); Radio bursts (1339); Magnetars (992) 


1. INTRODUCTION 
In 2007, Lorimer et al. (2007) reported FRB 010724 as the first detected fast radio burst (FRB). It was discovered 


during an analysis of archived pulsar survey data, recorded simultaneously by three adjacent beams on the receiver 
of Parkes radio telescope (Parkes) in Australia. It has a high Galactic latitude, but with a relatively low signal-noise 
ratio. The second FRB, FRB 010621, was discovered by Keane et al. (2012), which has a lower galactic latitude 
(b = —4.003?). Subsequently, Thornton et al. (2013) reported a sample of four new FRBs. Since then, the study of 


FRBs has become an important new field in high energy astrophysics. 
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The observation of FRBs has entered a new era since 2017. Nearly 10? FRB sources have been detected so far. 


However, due to the limited field of view of radio telescopes, the observed event rate is far less than the intrinsic FRB 
rate in the whole sky. FRBs can be divided into two main groups, i.e. repeating FRBs and one-off (non-repeating) 
FRBs. 5o far, the number of observed one-off FRB sources is much more than those of repeating FRBs. With the 
accumulation of observational data, it is getting more and more clear that FRBs should be connected with compact 
objects, especially highly magnetized neutron stars. Viable models include: giant flares from magnetars (Popov & 
Postnov 2010; Kulkarni et al. 2014a,b), giant pulses of young pulsars (Connor et al. 2016; Cordes & Wasserman 
2016; Pen & Connor 2015), star quakes of young magnetars (Lieu 2017), the interaction between magnetar flares 
and their surroundings (Lyubarsky 2014; Katz 2016a; Beloborodov 2017; Metzger et al. 2017b), lightning of pulsar 
magnetospheres (Katz 2017), planets around pulsars (Mottez & Zarka 2014), collapse of massive neutron stars to 
black holes (Falcke & Rezzolla 2014; Zhang 2018a; Ravi & Lasky 2014), neutron star-white dwarf/Be star binaries (Gu 
et al. 2016; Li et al. 2021c), binary neutron stars mergers (Totani 2013; Wang et al. 2016), collisions between neutron 


stars and asteroids (Geng & Huang 2015; Dai et al. 2016). Other novel models even include: cosmic string collisions 
(Cai et al. 2012), stellar flares (Loeb et al. 2014), binary white dwarf mergers (Kashiyama et al. 2013), evaporation 
of primary black holes (Barrau et al. 2014), mergers of charged binary black holes (Zhang 2016a; Liu et al. 2016; 
Punsly & Bini 2016), collapses of the crust of strange stars (Geng et al. 2021), light sails of extragalactic civilizations 
(Lingam & Loeb 2017), collisions of relativistic field structures (Thompson 2017), and interactions of strong plasma 
stream with pulsar magnetospheres (Zhang 2017). In this new field of high energy astrophysics, further exploration 
of FRBs will help us better understand binary neutron star mergers (Li et al. 2021c; Callister et al. 2016), supernova 
remnants (Yang et al. 2016; Piro 2016; Metzger et al. 2017a; Li et al. 2020a), cosmic baryon content (Macquart et al. 
2020a; Li et al. 2020b), compact faint objects (Munoz et al. 2016; Sammons et al. 2020), new cosmic ray acceleration 
mechanisms (Litwin & Rosner 2001). It will also help examining basic physics (Wei et al. 2015; Wu et al. 2016; Zhang 
2016b; Tingay & Kaplan 2016; Nusser 2016; Shao & Zhang 2017) and constraining cosmology models and dark energy 
content (Gao et al. 2014a; Zhou et al. 2014a; Yang & Zhang 2016; Yu & Wang 2016, 2017; Li et al. 2018; Zhao et al. 
2020). 

It is interesting to note that FRBs could be connected with some kinds of X-ray bursts (XRBs) from magnetar. In 
the history, magnetars were initially proposed to explain soft gamma repeaters (SGRs) by Duncan & Thompson (1992). 
They are the most extreme type of pulsars characterized by very strong magnetic fields, being born as the result of the 
evolution of massive stars towards the end of their lives or the binary neutron star mergers. Magnetars continuously 
emit high energy electromagnetic radiation, mainly in X-rays and ^-rays, through the decay of its magnetic field. 


Originating in magnetars, SGRs enter an active phase every few years or decades, generating a series of X-ray and 


y-ray bursts. Recent observations have revealed that at least some FRBs are related to XRBs from magnetars, but 
the detailed mechanisms are still unknown. 

To be specific, some repeating FRBs have periodic window behavior (PWB) that active periods occur quasi- 
periodically (CHIME/FRB Collaboration et al. 2020; Rajwade et al. 2020b; Cruces et al. 2021). As a valuable SGR 
sample found to be closely related to FRBs, SGR J1935+2154 was first detected by the Swift Burst Alert Telescope 
(BAT) in 2014 (Stamatikos et al. 2014). FRB 20200428A was observed to be associated with this active SGR by 
the Hard X-ray Modulation Telescope (HXMT) and the Canadian Hydrogen Intensity Mapping Experiment telescope 
(CHIME) (Younes et al. 2017; Lin et al. 2020), in association with one of a series of XRBs (CHIME/FRB Collaboration 
et al. 2020; Li et al. 2021a; Mereghetti et al. 2020; Ridnaia et al. 2021; Bochenek et al. 2020). Such an association 
explicitly supports the idea that at least some FRBs might be connected with magnetars. Inspired by FRB 202004284, 
researchers have extensively investigated the relation between periodically repeating FRBs and PWBs of SGRs (Zhang 
et al. 2021; Grossan 2021; Denissenya et al. 2021; Zou et al. 2021). 
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However, since FRB 20200428A is the only sample observed to be associated with an XRB so far, the origin and 


trigger mechanism of FRBs are still highly debated. In order to better reveal the nature of FRBs, here we present a 
thorough study on repeating FRBs. A large sample of repeating FRBs observed by 11 different telescopes are collected. 
The data are analyzed and compared for each FRB source, aiming to provide useful guidance for future observation 


strategies and modelling of repeating FRBs. 


2. DATA AND METHODS 


Till now, more than 700 FRB sources have been detected, 29 of which are repeating FRBs (see chime-frb-web, 
HeRTA: FRBSTATS, Transient Name Server, FRB Catalogue and Fast Radio Burst Community Newsletter). It is not 
clear whether different repeating FRBs have the same physical origin, but it is interesting to see whether the observed 
parameters of repeating FRBs are similar or not. We have collected the observational data of repeating FRBs from the 
literature, and will present a statistical analysis on them. Four FRB sources, i.e. FRB 20121102A, FRB 20180916B, 
FRB 20190520B and FRB 20201124A are the most active repeating events, which will be paid special attention on in 
this study. 

FRB 20121102A, as the first detected repeating FRB, has the longest observation history. It was first discovered 
by the 305-m William E. Gordon Telescope at the Arecibo Observatory (Arecibo) (Spitler et al. 2014) and was later 
followed by a large number of observations (Spitler et al. 2016; Hewitt et al. 2022b; Michilli et al. 2018; Hilmarsson 
et al. 2021b). More than 2000 FRBs have been detected from this source. Especially, between August 29 and October 
29, 2019, 1652 bursts were detected by the Five-hundred-meter Aperture Spherical radio Telescope (FAST) (Li et al. 


2021b). The average dispersion measure (DM) is 565 cm ? pc, and the average Faraday rotation measure (RM) is 


78153 rad m ?. The latest analysis suggests that there is no periodicity in its activity (Katz 2022). It located in 
an irregular dwarf galaxy with a low metallicity, at a redshift of z — 0.193. The source is in a star-forming region, 
associated with a luminous persistent radio source. The properties of the persistent radio source are consistent with 
a low luminosity, accreting massive black hole (Michilli et al. 2018; Chatterjee et al. 2017; Kokubo et al. 2017). 

FRB 20180916B was first discovered by CHIME (CHIME/FRB Collaboration et al. 2019), with about 90 bursts being 
observed initially. Between April and December 2020, 54 bursts were further observed by the Westerbork Synthesis 
Radio Telescope (WSRT), which consists of twelve 25-m dishes in a new system called the Aperture Tile in Focus 
(Apertif) (Pastor-Marazuela et al. 2021). The total number of bursts detected from FRB 20180916B exceeds 150. The 
average DM is 349 cm? pc and the average RM is -108 rad m~?. This FRB souce may have an active periodicity 


of 16.35 days (Chime/Frb Collaboration et al. 2020; Voisin et al. 2021). It located in a nearby massive spiral galaxy 
at a redshift of z — 0.0337. The source is in a star-forming region, but is not associated with any persistent radio 
counterpart (Marcote et al. 2020). 

FRB 20190520B was first discovered by FAST, with four bursts being reported initially. Subsequently, 75 bursts 
were further observed by FAST between April 25 and September 19, 2020 (Niu et al. 2022). Between June 2021 and 
January 2022, Parkes detected 113 additional bursts (Dai et al. 2022). The total number of bursts being observed 
from this source is close to 200. The average DM is 1207 cm"? pc, and the average RM is 3603 rad m~?. This 
FRB has a very distinct DM excess and no periodicity has been found in its activity (Niu et al. 2022). It located in 


a dwarf galaxy with a high specific star formation rate (SFR) at a redshift of z = 0.241. The source is in an extreme 
magneto-ionic and dense environment, in association with a compact persistent radio source. Such an environment 
may hint a distinctive origin or an earlier evolutionary stage for this FRB source (Niu et al. 2022). 

FRB 202011244, as a newly detected active repeating FRB, was first discovered by CHIME. CHIME detected more 
than 33 bursts from it (chime-frb-web). From April 1 to June 11, 2021, FAST has detected 1863 bursts from this 
source (Xu et al. 2022). The total number of bursts observed from FRB 20201124A thus exceeds 2000. The average 
DM is 413 cm ^? pc, and the average RM is -586 rad m^?. The RM of this FRB was found to vary significantly over 
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time (Wang et al. 2022a). No periodicity is found in its activity (Xu et al. 2022). It locates in a Milky-Way-sized 
barred spiral galaxy with a high metallicity, at a redshift of z = 0.09795. ‘The source resides in an inter-arm region at 


an intermediate distance from the center of its host, with a low stellar density. Such an environment does not favor 


a young magnetar origin, which requires an extreme explosion of a massive star (Xu et al. 2022; Nimmo et al. 2022c; 
Ravi et al. 2022). 

For these active repeating FRBs, statistical studies have been conducted extensively on individual FRBs. For 
instance, studies on the periodicity (Katz 2022), spectral features (Levkov et al. 2022; Rajwade et al. 2020a), radio 
variability (Chen et al. 2022a), and parameter statistics (Hewitt et al. 2022a; Li et al. 2019) are carried out for FRB 
20121102A. For FRB 20180916B, periodicity (Chime/Frb Collaboration et al. 2020; Voisin et al. 2021) and multi-band 
alignment (Sand et al. 2022) are investigated. Also, RM variation of FRB 20190520B (Dai et al. 2022; Anna-Thomas 
et al. 2022) has been explored, and polarization variation (Kumar et al. 2022) and parameter statistics (Marthi 
et al. 2022; Lanman et al. 2022) are conducted for FRB 20201124A. Meanwhile, some researchers have studied the 
classification of FRB populations through various features connected to pulse profile, brightness temperature (Chaikova 
et al. 2022), DM (Chawla et al. 2022), location within host (Chrimes et al. 2021), SFR (James et al. 2022), source 
density (Law et al. 2022), and parameter synthesis (Gardenier et al. 2019, 2021). Notably, there is a dominant view 


that repeating FRBs and one-off FRBs may be two distinct populations at all (Pleunis et al. 2021; Rafiei-Ravandi et al. 
2021; Zhong et al. 2022; Chen et al. 2022b). With the accumulation of observational data, new statistical researches on 
various parameters are updated (Katz 2016b; Macquart & Ekers 2018b,a; Ravi 2019; Hashimoto et al. 2020a; Takahashi 
et al. 2021; Nimmo et al. 2022a; Bhandari et al. 2022), constantly updating our understanding of these interesting 
sources. Note that there are also some statistical studies that concentrate specifically on one-off FRBs (Hashimoto 
et al. 2020b, 2022) or repeating FRBs (CHIME/FRB Collaboration et al. 2019; Feng et al. 20222). 


A recent study by Zhang et al. (2022a) suggests that there are some statistical similarities between one-off and 


repeating FRBs. However, note that it is actually a very difficult issue to compare repeating and one-off FRBs 
directly. Repeating sources are usually constantly monitored at high sensitivity, which can lead to the detection of a 
large number of weak FRBs. As a result, the distribution of parameters may be systematically different for repeating 
FRBs and one-off FRBs. In this aspect, it has been suggested that the difference between one-off FRBs and the first 


detected burst of repeating sources should be somewhat comparable and could present meaningful information on their 


origin. Additionally, it is still controversial as to whether the one-off FRBs actually repeat or not. In response to these 
issues, we will limit our research to repeating FRBs only. We have collected the observational data from 11 different 
telescopes, including FAST of China, CHIME in Canada, Arecibo in USA, the 76-m Lovell telescope (Lovell) in UK, 
the Effelsberg 100-m Radio Telescope (Effelsberg) in Germany, the Apertif in Netherlands, the LOw Frequency ARrray 
(LOFAR) in Europe, the upgraded Giant Metrewave Radio Telescope (uGMRT) in India, the Parkes in Australia, the 
Robert C. Byrd Green Bank Telescope (GBT) and the KarlG. Jansky Very Large Array (VLA) in USA. The data 
correspond to various repeating FRB sources, including multiple frequency observations at different times. An overall 
view of the data sources are shown in Table 1. 

Various statistical analyses are conducted on these repeating sources. Unbiased estimation under normal or lognormal 
distribution is used to standardize the statistical characteristics of each parameter for each repeater. Single-parameter 


analyses (Section 3) and double-parameter analyses (Section 4) are carried out to explore intrinsic features of repeating 


FRBs. Based on these investigations, conclusions are presented for each source (Section 5). 


3. SINGLE-PARAMETER ANALYSIS 


The emission of FRBs is mainly in radio wavebands (GHz), lasting typically for a few milliseconds with the flux 


density being of the order of several Jansky. FRBs have significant extragalactic DM excess. Such an excess is caused 
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Table 1. An overview of the main data used in this study. The data were collected as of June 26, 2022. In the top half 
of the table, the number of bursts detected by different telescopes is presented for each repeating FRB source. Note that 
FRBs 20121102A, 20180916B, 20190520B and FRB 20201124A are the most active sources. There are still some other inactive 
repeating FRBs with less than 10 bursts being observed (mainly observed by CHIME). They are not listed in this table but 


are used for DM analysis below. The bottom half of the table indicates the parameters available from each telescope. 


FRBname FAST CHIME* Arecibo  Lovell Effelsberg Apertif LOFAR uGMRT Parkes GBT* VLA "Total 
20121102A 1652 518 32 4 93 2 2301 
20201124A 1863 34 20 48 1965 
20190520B 79 113 3 195 
20180916B 91 54 9 154 
20200120E 8 60 68 
20190303A 27 27 
20180814A 22 22 
20180301A 20 1 21 
20181119A 13 13 
20190417A 13 13 
20190208A 11 11 
20190212A 11 11 

MJD* Vo xw ow Vv v v v v VEA d v v u^ y v V o s 

DM Vo v xv X v v d v y y VA v VA 
RM V v/ Part V v/ vA V Lo EA 

Width V v v x Part V V V vw sy V A vA 
Bandwidth v "d Part d v v v "d 
Peak Flux vV d 74 Part d v v v NA v A 

Fluence Y VY vw ^ Part v v v EVA v v d v v v 


Reference (1) (2,3) (4) (5) (6,7 (8,9,10,11) (12) (11,13) (14) (15) (16) (16) (17) (18) (19) (20,21) (22) (11) 


NOTE— References: (1) Li et al. (2021b); (2) Luo et al. (2020); (3) Laha et al. (2022); (4) Niu et al. (2022); (5) Xu et al. (2022); (6) CHIME/FRB 
Collaboration et al. (2021); (7) Mckinven et al. (2022); (8) Spitler et al. (2016); (9) Hewitt et al. (2022b); (10) Michilli et al. (2018); (11) 
Hilmarsson et al. (2021b); (12) Rajwade et al. (2020b); (13) Spitler et al. (2018); (14) Nimmo et al. (2022b); (15) Hilmarsson et al. (2021a); (16) 
Pastor-Marazuela et al. (2021); (17) Marthi et al. (2022); (18) Price et al. (2019); (19) Dai et al. (2022); (20) Gajjar et al. (2018); (21) Zhang 
et al. (2018); (22) Feng et al. (2022b). 

* CHIME provides MJD and DM data for each burst, but many other physical parameters are available only for part of these bursts. 
+ FRB 20121102A was observed by GBT, but note that the parameters reported in these two references are slightly different. 
* MJD is useful in calculating the waiting time. 


by the intergalactic medium (IGM), suggesting that FRBs should be of extragalactic or even cosmological origin 
(Deng & Zhang 2014a; Gao et al. 2014b; Kulkarni et al. 2014c; Zheng et al. 2014; Zhou et al. 2014b). Direct spectral 
observations of the optical counterpart of FRB 20121102A show that its redshift is z = 0.19273 + 0.00008, confirming 
its extragalactic origin(Tendulkar et al. 2017). These facts mean that a typical FRB can release an isotropic energy up 
to 10*! erg. The short duration of FRBs implies that the source should be small and may be connected with compact 
stars. The intensity distribution function can also provide useful clues on the origin of FRBs (Li et al. 2017). It has 
also been noted that the galactic latitudes of FRBs are usually high (ranging from 20° to 70°). By analyzing the dipole 
moments (cos @) and quadrupole moments (sin? b — 1/3) of FRB’s distribution on the sky, it is found that FRBs are 
quite isotropic (Keane & SUPERB Collaboration 2016). 

Considering the intense emission in radio frequencies, the brightness temperatures of FRBs are generally as high as 
~ 10?? K, i.e. 


F sa pm ure 
2 2 35 V 
Tg ~ F, ,D?/2rk (vT)? z 1.1 x 1095 K (3 (= -) (Z) (3 (1) 
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where, I, = F,/n10? is the radiation intensity at frequency v, 0 represents the angular radius of the source, k is the 
Boltzmann constant, and A is the wavelength. 0 can be approximated as cT'/ D, where c is the speed of light, T' is 
the time scale of the FRB, and D is the distance of the source. The high brightness temperature of FRBs indicates 
that non-thermal radiation should be involved (Luan & Goldreich 2014; Katz 2014). The light curves of FRBs mainly 
show a single pulse structure, but Champion et al. (2016) found a bimodal pulse profile in one burst from the source 
of FRB 201210024. Recently, Pastor-Marazuela et al. (2022) also reported that FRB 20201020A has a sub-millisecond 


quasi-periodic structure. FRBs are mainly characterized by several key parameters, such as DM, RM, bandwidth, 


pulse width, waiting time, peak flux and fluence. In this section, we present a detailed analysis on these parameters. 


3.1. Dispersion Measure 


When the electromagnetic wave propagates in plasma, it interacts with free electrons and becomes dispersed. The 


group velocity of electromagnetic waves becomes frequency-dependent. As a result, high frequency emission arrives 


earlier than low frequency waves. The dispersion measure is defined as the integral of electron number density (ne) 


along the line of sight, which can be expressed as 


DM = / yay = (Ne) D (2) 


Obviously, DM can be divided into several parts (Xiao et al. 2021), 


DM Hos 
DM = DMuw + DMiem + [—. (3) 


Here, DM mw is the contribution from the interstellar medium (ISM) within the Milky Way (MW), DM ray is the 
contribution of the intergalactic medium (IGM), and DM Host corresponds to the contribution of ISM within the host 
galaxy and the local circum-burst environmental near the source. 
The arrival time of the electromagnetic waves is frequency dependent, and the delay time at frequency v is related 
to DM as (Wilson et al. 2009), 
At, ee i UCL (4) 


QTMeC y? 


where e is the electron charge, and m, is the electron mass. Using this equation, we can estimate the DM of FRBs 


through multi-wavelength radio observations, 


29 
2TMeC UrLU| 


2 (5) 


DM = (Ati — Ath) a a 
h l 


where Af; represents the time delay of the pulse at a lower frequency nı, and At, represents the time delay of the 


pulse at a higher frequency vp. Sometimes, multipath propagation caused by scattering makes the plasma act like a 
lens, and the DM can be used to establish plasma lens models. 

DM is an important parameter of FRBs and is usually used as a distance probe. Firstly, we can estimate DM yw 
by using the widely used MW free electron models such as NE2001 (Cordes & Lazio 2002) or YMW16 (Yao et al. 
2017). Through the electron distribution after the Epoch of Reionization (Deng & Zhang 2014b), we can also calculate 
DM em at redshift z (Inoue 2004; Ioka 2003). Combining DM Host estimated by assuming a particular galaxy type 
(Luo et al. 2018), we can finally infer the redshift and luminosity distance of the source by Formula (3). Moreover, 
taking into account the visibility of various observing epochs can also help to estimate the distance of the FRB source 
(Nimmo et al. 2022c). On the other hand, it is also possible for us to estimate DM ram through the Macquart relation 
(i.e. DM raw — z relation) when the redshift is available for the repeating FRB (Macquart et al. 2020b). 
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Figure 1 plots the DM of FRBs as a function of time. We see that in most cases, DM does not have an obvious 
evolution trend over time. However, there are still several FRBs seeming to be special. The DM of FRB 20190520B 
shows a rapid declining trend, while the DM of FRB 20201124A and FRB 20121102A show an increasing trend (also 


see the insets in Figure 1). For the four most active FRBs, we have fitted their DM evolution by using a linear function, 


and then pay special attention on the DM residuals with respect to the linear function. The residual (dispersion of 
DM ) reflects the random variation range of DM. The residuals are also plotted in Figure 1. We see that the DM 
residuals generally follow a normal distribution for all the four sources. It is interesting to note that the DM variation 
is mainly in a range of +3 cm ? pc. However, FRB 20190520B seems to be special since its variation is as large as 
+12 cm ^? pc. It is still unclear whether such a variation is due to the change of the local circum-burst environment 


Or not. 


In Figure 2, various correlations between redshift, SFR and DM are plotted for repeating FRBs. In Panel (a), we 
see that DM is linearly correlated with redshift. It is consistent with the results of Macquart et al. (2020b), which 
reflects the major role played by DM jay in DM. The root mean square error (RMSE) is used to characterize the 
dispersion level of DM with respect to the linear fit for each repeating FRB. In Figure 2(b), we see that DM RMSE is 
also positively correlated with redshift, which may be due to IGM fluctuations. In Figure 2(c), DM increases linearly 
with the SFR of the host galaxy. It should be caused by the increase of DM Host thanks to star formation activities 
in the host galaxy. In Figure 2(d), the RMSE and time-varying rate of DM are both normally distributed. There is 
also no correlation between these two parameters. Figure 2(e) shows that the time-varying rate of DM is not related 
to redshift. In all the five panels of Figure 2, we notice that FRB 20190520B is very different. It is always an outlier 
with respect to other FRBs. This could be caused by extreme activities in the source or host galaxy, or by strong 
turbulence on the line of sight (Anna-Thomas et al. 2022). 

The DM of FRB 20190520B decreases rapidly over time, while the DM of FRB 202011244 increases rapidly. The 
underlying physical mechanism could be quite different for them. Note that the DM abnormality in FRB 20190520B 


was initially reported by Niu et al. (2022). Here, we have carried out further statistical analysis on this issue. Detailed 


statistical characteristics are presented in Table 2. 


3.2. Faraday Rotation Measure 


RM is another key parameter of FRBs and is often used to infer the magnetic field on the line of sight. Due to the 


Faraday effect, the polarization angle Woss (A) at wavelength A can be expressed as (Xiao et al. 2021) 


V ,5.(A) = Voc AW = Vo + 


e3A? [ Ne (1) Bi (1) 7 (6) 


iwi Jo [xig 


where Wo is the initial polarization angle, AW is the variation of the polarization angle induced by the Faraday effect, 


and Bj represents the magnetic field parallel to the line of sight. RM is defined as 


3 P n, (1) By (I 
RM — aa | ne BV) (7) 
Irme Jo ez) 
Then Yobs (A) can be simplified as: 
U 
Voss (A) = Vo +A° RM = 0.5 arctan Q (8) 


obs (A) can be calculated from the observed Stokes parameters U and Q, so that RM can also be derived through 
observations (Hotan et al. 2005; Macquart et al. 2012; O'Sullivan et al. 2012). Furthermore, combining DM and RM, 
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Figure 1. DM plotted versus time for each repeating FRB source. The horizontal dashed line corresponds to the DM average 
of each FRB, and the solid line is the linear fit to the observational data points. The inset shows the observational data points 
from different telescopes. For the four most active FRBs, the residuals of data points with respect to the linear fit are also 


illustrated via histogram. 


we can estimate the parallel magnetic field as (Caleb et al. 2018) 


RM RM DM X! 
(B1) = 0.812DM ~ ads rad m~? pe em-? i (9) 


Note that the linear polarization of FRBs will be significantly affected during the propagation. A recent study by Feng 


et al. (2022a) showed that FRBs depolarize at lower frequencies, i.e. the degree of linear polarization decreases with 


frequency. 
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Figure 2. Correlations between redshift, SFR and DM for the repeating FRBs. Panel (a): the relation between DM and 
redshift. Panel (b): the relation between DM RMSE and redshift. Panel (c): the relation between DM and SFR. Panel (d): 
the relation between DM RMSE and DM time-varying rate. Panel (e): the relation between DM time-varying rate and redshift. 


Figure 3 plots RM versus time for all the repeating FRBs with RM data available. We see that RM generally has a 
strong evolution over time (Dai et al. 2022; Wang et al. 2022a). This can be seen in all the five FRBs shown in Figure 
3, i.e. FRBs 20121102A, 20180301A, 20190520B, 202011244, and 20201124A. However, the evolution trend is different 
for different FRBs. In some FRBs, RM increases with time, but in others, RM decreases with time. In the case of 
FRB 20201124A, which has the most prolific RM data, the evolution is very complicated. We notice that the RM of 
FRB 20180301A and FRB 20190520B experienced variations between positive and negative values. It could be caused 


by the inversion of the magnetic field in the emission region, but it is also possible that FRBs could happen in more 


than one region near the source. To some extent, the former case does not support the young magnetar model unless 
there is a special magnetization mechanism; while the latter may reveal the involvement of orbital movements. For 
example, the triggering position may vary between the equator and the polar regions. Whether there is any periodicity 
in the variation of RM still needs to be examined. In the histogram plots, we can clearly see the variation range of 
RM in each repeating source. In the case of FRB 202011244, the histogram can be fitted by a normal distribution, 
but with a high peak at the center. 


The mean value of RM can be regarded as an indication of the magnetic field strength. Here we use |RM| to refer 
to its absolute value. Figure 4 shows the relations between redshift, SFR and | RM| for repeating FRBs. We see that 
|RM| has no correlation with redshift, but it correlates negatively with SFR. This could support young magnetars 


as the central engines of repeating FRBs. Combining with the above timing evolution of RM, we argue that the 


triggering position of bursts should be constantly changing near the source. Detailed statistical analyses are shown in 


Table 3. 
Figure 5 shows the correlation between the statistical parameters of RM for repeating FRBs. We see that the 


statistical parameters, i.e. the average level, dispersion level and peak of RM distribution, are positively correlated 
with each other. However, it should be noted that FRB 20190520B seems to be an outlier in these plots. This is 
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Table 2. Statistical characteristics of DM for all FRBs with an observed host galaxy. 
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DMZ 


Type FRBname Redshift References” DM DM Rate DM RMSE Exc DM... SFR RM 
(cm? pc) (cm ? pecyr !) (cm ? pc) (cm ? pc) (cm ? pc) (Mo yr !) (radm ?) 
Repeating 20121102A 0.1930 (1,2,3) 565.83 1.76 1.95 347.80 180.17 78152.95 

20180916B 0.0337 (4) 349.20 -0.26 2.44 150.40 121.13 -107.83 
20190520B 0.2410 (5) 1207.04 -7.55 10.68 1097.04 887.72 0.41 3603.18 
20201124A 0.0980 (6,7,8) 413.27 5.67 2.14 273.27 188.20 1.70 -586.59 
20180301A 0.3304 (9) 519.44 -2.73 3.52 367.44 80.48 1.93 81.96 
20181030A 0.0500 (10) 103.70 0.12 1.93 62.30 18.87 0.36 
20200120E 0.0008 (11,12) 88.20 -0.95 0.24 48.20 47.51 0.40 - 0.80 
20190711A 0.5220 (13) 593.10 

One-off 20150418A 0.4920 (14) 776.20 
20180924B 0.3212 (15) 361.42 
20181112A 0.4755 (16) 589.27 
20190102C 0.2912 (9) 363.60 
20190523A 0.6600 (17) 760.80 
20190608B 0.1178 (18) 338.70 
20190611B 0.3778 (9) 321.40 
20190714A 0.2365 (9) 504.00 
20191001A 0.2340 (19) 506.92 
20191228A 0.2430 (9) 297.90 
20200430A 0.1600 (9) 380.10 
20200906A 0.3688 (9) 577.80 


NOTE— References: (1) Michilli et al. (2018); (2) Chatterjee et al. (2017); (3) Kokubo et al. (2017); (4) Marcote et al. (2020); (5) Niu et al. (2022); (6) Xu 
et al. (2022); (7) Nimmo et al. (2022c); (8) Ravi et al. (2022); (9) Bhandari et al. (2022); (10) Bhardwaj et al. (2021a); (11) Bhardwaj et al. (2021b); (12) 
Kirsten et al. (2022); (13) Heintz et al. (2020); (14) Keane et al. (2016); (15) Bannister et al. (2019); (16) Prochaska et al. (2019); (17) Ravi et al. (2019); 
(18) Chittidi et al. (2021); (19) Bhandari et al. (2020). 

* Redshift and SFR data were obtained from these references. 
+ The extragalactic DM is calculated by using the NE2001 model (Cordes & Lazio 2002). 
* DM component contributed by the host galaxy is calculated by using the Macquart relation (Macquart et al. 2020b). 


mainly due to its large |RM| value and the DM abnormality possibly produced by special activities in the source or 


its host galaxy, as discussed in Section (3.1). 


3.3. Bandwidth 


FRB emission is usually in a narrow frequency range in the radio band. The observed bandwidth is an important 
parameter of FRBs, which is also limited by the work band of telescopes. For this reason, we need to analyze the 
observational data from each telescope separately. 

Figure 6 plots the bandwidth versus time (MJD) for the most active repeating FRBs. The distribution of the 
bandwidth is also illustrated via histogram. We see that the bandwidth does not show any obvious time evolution 
trend and it is generally normally distributed. For the Arecibo data of FRB 20121102A, it can be seen that fj; and 
frow stack up at the upper and lower bounds of Arecibo’s work band, respectively (see Figure 6, Row 1, Column 3). 
It is due to the cutoff of the limiting work band of Arecibo. The intrinsic bandwidths of these bursts should be beyond 
the work band. Similar phenomenon can also be seen in other panels, such as in Row 1, Columns 4 and 3). Note that 


the limiting work band of telescopes could also lead to the distortion in the fluence data (Keane & Petroff 2015). 
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3. Evolution of RM for repeating FRBs. The horizontal dashed line indicates the average value of RM for each FRB. 


The solid line is the best linear fit to the data points. For each FRB, the distribution of RM is illustrated via histogram. The 


inset shows the observed data points of each telescopes. 
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Figure 4. Correlations between redshift, SFR and RM for repeating FRBs. Panel (a): the relation between | RM | and redshift. 
Panel (b): the relation between |RM| and SFR. Note that the linear fit in this panel very tight so that the 90% confidence 


interval is too narrow to be seen. 
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Figure 5. Correlation between the statistical parameters of RM for repeating FRBs. u and o represent the two parameters 


of normal distribution of RM for each FRB, characterizing the average value and the dispersion level, respectively. Average 


stands for the mean value, which is different from p here because p is derived by a normal distribution fitting. Peak represents 


the value observed most frequently. 


Table 3. Detailed statistical characteristics of RM for repeating 
FRBs. The table contains the main statistical parameters about the 
normal distribution of RM for each FRB. u and ø represent the two 


parameters of normal distribution, characterizing the average value 


and the dispersion level, respectively. Average stands for the mean 
value, which is different from u here because u is derived by a nor- 


mal distribution fitting. Peak represents the value observed most 


frequently. 

FRBname u c Average Peak 
(rad m7?) (rad m7?) (rad m7?) (rad m7”) 

20121102A 69204.60 12365.60 18152.95 69204.63 

20180301A 81.96 1226.54 81.96 191.21 

20180916B -112.90 20.50 -107.83 -112.91 

20190520B 9641.80 10016.06 3603.18 9641.80 

20201124A -586.59 105.72 -586.59 -601.85 


Figure 7 plots the characteristic bandwidth versus redshift for repeating FRBs. Here the parameter UBandwidth 
can be regarded as the typical bandwidth of an FRB source. We see that Upandwiath has a clear positive correlation 
with redshift. Note that Arecibo (1.15-1.73 GHz), FAST (1.05-1.45 GHz) and Effelsberg (1.2-1.6 GHz) have similar 


work bands. Their data are employed in Figure 7 to minimize the selection effect caused by instrumentation. Such a 
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correlation may be induced by the cosmological broadening of the bandwidth, but it could also be due to cosmological 
evolution of the central engines. 

Detailed statistical parameters of the bandwidth are listed in Table 4. Figure 8 shows the correlation between these 
parameters for repeating FRBs. We see that the parameters of bandwidth are positively correlated with each other 
except for the Peak-o plot. Note that these correlations are generally very weak (except for the -Average plot). A 


larger sample size will help clarify the correlations in the future. 
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Figure 6. Time evolution of the bandwidth for the active repeating FRBs. The horizontal dashed line corresponds to the 
average of the bandwidth for each FRB. The solid line is our linear fit to the data points. For the FRBs detected by different 
telescopes, the distribution of the bandwidth is illustrated via histogram. fHigh and fLow stand for the upper and lower frequency 


of the observed bandwidth, respectively. 


3.4. Pulse Width 


Pulse width, or the duration of a burst, is strongly influenced by the plasma on the line of sight when an FRB 
propagates toward us. Multipath propagation caused by scattering broadens the pulse and attenuates the peak flux. 
It also leads to a delay in the arrival time. The broadened pulse width of the FRB varies notably with the frequency, 
leading to a frequency-dependent pulse trailing (Williamson 1972; Lee & Jokipii 1975; Rickett 1977) as 


wx DM?-v^*. (10) 
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Figure 7. Correlation between redshift and the characteristic bandwidth for repeating FRBs. Bandwidth is the characteristic 


bandwidth derived from the normal distribution fitting, which could represent the average bandwidth. 
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Figure 8. Correlation between the statistical parameters of bandwidth for repeating FRBs. u and o represent the two 
parameters derived from the normal distribution fitting of the bandwidth for each source, characterizing the typical bandwidth 


and the dispersion range, respectively. Average stands for the mean value. Peak represents the value observed most frequently. 


The measured width is affected by the sensitivity and threshold of telescopes, so we analyze the observational data of 
each telescope separately. 

Figure 9 plots the pulse width versus MJD time for repeating FRBs. We see that the width does not have an apparent 
variation over time. It generally satisfies the lognormal distribution well. Lg, is derived in such a lognormal fit, 
which represents the typical pulse width for the FRB source. Detailed parameters derived are listed in Table 5. Figure 


10 illustrates the correlation between redshift, SFR and pulse width for repeating FRBs. There is a positive correlation 
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Table 4. Detailed characteristics of the bandwidth for repeating 


FRBs. u and o represent the two parameters derived from the nor- 
mal distribution fitting of the observed bandwidth for each source, 
characterizing the typical bandwidth and the dispersion range, re- 
spectively. Average stands for the mean value. Peak represents 


the value observed most frequently. 


FRB " c Average Peak 

(MHz) (MHz) (MHz) (MHz) 

20121102A Arecibo 295.26 122.04 295.26 265.44 
20121102A FAST 200.73 112.19 200.73 104.11 
20200120E Effelsberg 180.65 38.84 180.65 188.85 
20201124A FAST 211.26 111.15 211.26 167.21 


between Uwiath and redshift. It may be caused by a combination of the cosmological expansion and the propagation 
effects. Cosmological expansion will broaden the width at high redshift, and the scattering in a longer path will also 
broaden the width. However, note that cosmological evolution of FRB sources may also play a role in the correlation, 
which requires further investigation. FRB 20201124A observed by uGMRT seem to be an outlier in the uwidth — z 
diagram, the reason of which still needs to be clarified. Besides, there seems to be a weak positive correlation between 
LLwidth and SFR, but the data points are too few to draw any firm conclusion. Future observations will help to clarify 
it: 

Figure 11 illustrates the correlation between various statistical parameters derived by fitting the pulse width dis- 
tribution with a lognormal function for each repeating source. The parameters generally correlate with each other. 
Especially, the dispersion range of width is negatively correlated with the typical width itself, the dispersion range is 
negatively correlated with the peak, and the typical width is positively correlated with the peak. Note that the typical 
width and peak width of FRB 20201124A observed by uGMRT are higher than the fitting results. This is primarily 
caused by the higher detection threshold of fluence for uGMRT (Marthi et al. 2022). 


3.9. Waiting Time 


Waiting time is the time interval between two consequent bursts. It is closely related to the trigger mechanism of 
repeating FRBs, and is also an important factor that should be considered when modeling repeating FRBs. Note that 
the waiting time can be extracted only when the source is continuously monitored, which means the upper limit of the 
observed waiting time is subject to the observation duration. The wide field of CHIME makes it effective in measuring 
waiting times, since it has a relatively long monitoring period on any particular source, i.e., 

Üg-w 
360° 
where 0g. .w is the width of the east-west field of view of the telescope (in degrees), and 6 is the declination of the 
FRB source. For CHIME, 0g-w ~ 1.3? — 2.5°. According to Equation (11), the monitoring period of CHIME for 


FRBs near the north celestial pole can even be limitless. 


tobsery = 8616A4s - 


- Sec Ò, (11) 


Figure 12 plots the waiting time versus the MJD time for the most active repeating FRBs. We see that the waiting 
time does not have any evolution trend with time. The distribution of the waiting time is also illustrated via histogram. 
The sample size of the bursts detected by CHIME, Apertif and uGMRT is generally small (Marthi et al. 2022; Pastor- 


Marazuela et al. 2021). For all other samples, we could clearly see an interesting bimodal distribution in the histogram 


16 HU ET AL. 


$ FAST $ Arecibo $ WSRTApertif $  LOFAR $  Parkes = VLA == Average 
$  CHIME $ Lovell $  Effelsberg $ | uGMRT © GBT — Linear Fit 


FRB20121102A FRB20121102A Arecibo FRB20121102A FAST FRB20121102A Lovell 


= x T r 7 T 
=== Log-Norm Dist (u=1.4 o=0.6) F ! === Log-Norm Dist (u=1.4 o=0.6) 


=== Log-Norm Dist (u=0.9 o=0.5) | 
Average 


E 


== Average 


—- Average 


Counts 
Counts 


56500 57000 57500 58000 58500 109 10i LL L 4 n 
MJD Width (ms) Width (ms) aso 
FRB20LB0301A FRB20180301A FAST R Boa 
i ; i i * 10! FE Log-Norm Dist (u=1.6 070.4) r : J RE D T 7 T T — ——À : : i 
werage 


== Averag 


Counts 


Width (ms) 
Width (ms) 
Width (ms) 


= 


i 10° t l = 
L 1 L 1 L 1 i i i i i 
98200 58400 38600 58800 59000 59200 10! 58380 58400 58420 58440 58460 58480 56500 58600 58700 58800 58900 59000 59100 
MJD Width (ms) MJD MD 
FRB20190520B FAST FRB20200120E FRB20200120E Effelsberg FRB20201124A 
z T T T7 Y T T T T T T T T T T T T T 
m Log-Norm Dist (u=2.5 90.4) = == Log-Norm Dist (u=-1.7 0=0.6) ! - : 
—-— Average [ | == Average d 
* 
l + 
l B 
I 
10! | = à 10L 3 a 
g £ 2 1 E iL 
c ] = iL = 2 i = 
5 s 10 Í = i c 
3 ] $5 f a 9 
O = 9 ! ks] 
= fæ Z ! z 
L 
F I 
1 
1 
I 
i 
* 1 
10 I 
| ] : : E 1 1 1 1 1 1 19r m i i 1 1 1 1 1 , 
10! 59600 59610 59620 59630 59640 59650 3071 59310 59320 59330 59340 59350 59360 
Width (ms) MJD Width (ms) MJD 
FRB20201124A FAST FRB20201124A uGMRT 
=== Log-Norm Dist (u=2.0 0-04) = -- i i H = Log-Norm Dist (u=2.8 0=0.5) 
Average E he, 3 == Average 
LA aN 
[7 A N 
|) \ 10! 
"4 N 
N| u 
g 
4 =] 
o 
N à 
N 
K 
E 1 
10! 10? 
Width (ms) Width (ms) 


Figure 9. Pulse width plotted versus time for repeating FRBs. The horizontal dashed line corresponds to the average of the 
pulse width for each source. The solid line is our linear fit to the data points. The distribution of the data points is also 


illustrated via histogram. The inset shows the observational data points from each telescope. 


Champion et al. (2016), typically peaking at ~ 10 ms and ~ 100 s, respectively. Note that the events separated by 
~ 10 ms were argued by some authors as one single burst with a multiple-pulse sub-structure (Li et al. 2021b; Gourdji 
et al. 2019; Pastor-Marazuela et al. 2022; Chime/Frb Collaboration et al. 2022). Some of them may even be due to 
the strong gravitational lensing effect (Oguri 2019), or caused by the scintillation of IGM (Ravi et al. 2016) or ISM 
turbulence (Main et al. 2022). On the other hand, the waiting time of ~ 100 s should be close connected to the trigger 
mechanism of repeating FRBs (Lu et al. 2020). 

The parameter Lwaitingrime(High—value) 18 à measure of the typical waiting time of the right component in the 
bimodal distribution in Figure 12. It is derived by using a lognormal function to fit the right component for each source. 
The detailed results of such fits are presented in Table 6. Figure 13 illustrates the relation between redshift, SFR and the 


typical waiting time for repeating FRB sources. There is an obvious linear correlation between LwaitingTime(High—value) 


and redshift. Again, it may be induced by the joint effects of cosmological expansion as well as the cosmological source 
evolution. Note that FRB 20180916B is significantly higher than the fitting curve, which may be caused by the 
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Figure 10. Correlation between redshift, SFR and pulse width for repeating FRBs. Panel (a): Hwiatn versus redshift. Panel 
(b): Hwiath versus SFR. Note that uwiatn is derived by using a lognormal function to fit the observed pulse width distribution 
for each source, and thus represents the typical pulse width. FRB 20180916B is not used for statistical analysis due to the small 


amount of bursts, but is also shown in the figure for completeness. 
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Figure 11. Correlation between various statistical parameters derived by fitting the pulse width distribution with a lognormal 


function for repeating FRBs. u and o characterize the typical width and the dispersion range, respectively. Average stands for 
the mean value of pulse width. Peak represents the value observed most frequently. FRB 20180916B is not used for statistical 


analysis due to the small amount of bursts, but is also shown in the figure for completeness. 


relatively long duration of Apertif and CHIME observations (CHIME/FRB Collaboration et al. 2021). There is no 
correlation between Lwaitingrime(High—value) and SFR in Figure 13, which might imply that FRB arises from an 


external mechanism rather than an internal one (Li et al. 2019). 
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Table 5. Parameters derived by fitting the pulse width distribution 


with a lognormal function for repeating FRBs. u and o characterize 


the typical width and the dispersion range, respectively. Average 
stands for the mean value of pulse width. Peak represents the value 


observed most frequently. 


FRB u c Average Peak 
(In(ms)) (In(ms)) (ms) (ms) 

20121102A Arecibo 1.36 0.61 4.68 4.12 
20121102A FAST 1.38 0.64 4.90 3.21 
20121102A Lovell 0.87 0.47 2.66 2.53 
20180301A FAST 1.61 0.43 5.49 4.57 
20180916B CHIME* 0.70 0.66 2.50 1.06 
20190520B FAST 2.52 0.41 13.54 13.77 
20200120E Effelsberg -1.74 0.57 0.21 0.17 
20201124A FAST 1.95 0.45 7.77 T.11 
20201124A uGMRT 2.84 0.49 19.34 15.74 


NOTE— * With a small sample size. 


The correlation between various parameters concerning the waiting time are shown in Figure 14. The dispersion 


range of the waiting time is negatively correlated with the average value, and with the peak as well. The average value 
is positively correlated with the peak. In Figure 14 (Row 1, Column 1), the correlation between u and ø is different for 
the high-value and low-value waiting times, which indicates that the causes of them are different. Similar phenomena 
can also be seen in other two panels of (Row 2, Column 1) and (Row 2, Column 2). It is worth noting that FRB 
20190520B again is quite different from other sources in Figure 14. 


3.6. Peak Flux 


The peak flux, which represents the intensity of a burst, is widely used to calculate the fluence by multiplying it 
with the pulse width. It is a basic parameter for FRBs. 

The peak flux is plotted versus the MJD time for repeating FRBs in Figure 15. The distribution of the peak flux is 
also illustrated via histogram for each source. It could be seen that the peak flux has a weak decreasing trend on a 
long-time scale, which has been previously noted by Li et al. (2021b) on the bimodal burst energy distribution of FRB 
20121102A. In other words, strong bursts generally appear in an early stage, and late bursts are weaker in principle. 
Interestingly, for FRB 20121102A, Figure 15 (Row 1, Column 2) shows that the number of weak bursts is slightly 
larger than the fitting curve, which is different from that of FRB 20201124A (Row 2, Column 2, and Row 2, column 


3). In the future, when more FRBs are detected from these two sources, the histogram should be re-examined again 


to see whether it is still different for these two sources. 


In order to compare the characteristics of peak flux in different repeating sources, we still use a lognormal function 


to fit its distribution. The derived parameters are presented in Table 7. Here, UpeakFiur is a measure of the typical 
peak flux. UpeakFlux is plotted versus the redshift for repeating FRBs in Figure 16. It is still difficult to determine 
whether there is a correlation between uL pear riu; and redshift or not, due to the very limited data points in the figure. 


Correlations between various statistical parameters concerning the peak flux are shown in Figure 17. We see that the 


dispersion range of the peak flux is negatively correlated with the average peak flux. 
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Figure 12. The waiting time plotted versus MJD time for the most active repeating FRBs. The horizontal dashed line 
corresponds to the average value of waiting time for each source. The solid line is the best linear fit to the data points. The 
distribution of the waiting time is also illustrated via histogram. The inset shows the observational data points from each 


telescope. 


3.7. Fluence 


For an FRB source, the fluence represents the energetics of a burst. As mentioned earlier in Section (3.3), the 
measured fluence are affected by the work bandwidth, the detection threshold, and the instrument sensitivity of 
telescopes. The energy release of an FRB can be calculated through their fluence as 


4 
Ee J ^ — D? F,dv, (12) 


1+2z 


Where Dz is the luminosity distance, and F, is the fluence at frequency v. Zhang (2018b) suggested to use the center 
frequency of telescope work band to calculate the burst energy. Using this method, Li et al. (2021b) argued that the 
burst energy of FRB 20121102A satisfies a bimodal distribution. 

Figure 18 plots the fluence versus MJD time for repeating FRB sources. We see that the fluence generally does not 


have any obvious evolution on a long-time scale, which is slightly different to the behavior of peak flux. Comparing 
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Panel (a): 


LW aitingTime(High—value) V€rSus redshift. Panel (b): LWW aitingTime(High—value) VETSUS SFR. Here LW aitingTime(High—value) refers 
to the typical waiting time of the right peak in the bimodal distribution, derived by using a lognormal function to fit the peak 
for each FRB. FRB 201800301A and FRB 20181030A are not used for statistical analysis due to the small amount of bursts, 


but are shown in the figure for completeness. 


Table 6. Parameters derived by fitting the distribution of the waiting time with a lognormal function for repeating 


FRBs. u and c characterize the typical waiting time and the dispersion range, respectively. Average stands for 


the mean value of waiting time. Peak represents the value observed most frequently. 


FRB High — Value Low — Value 
" c Average Peak u c Average Peak 
(In(s)) (In(s)) (s) (s) (In(s)) (In(s)) (s) (s) 
20121102A Arecibo 4.44 1.01 141.54 101.89 -3.69 0.66 0.03 0.02 
20121102A FAST 4.11 1.09 111.02 72.20 -5.22 1.19 0.01 0.00 
20121102A Lovell* 6.09 0.82 617.50 320.28 - - - - 
20180301A FAST* 7.40 0.66 2037.87 1384.06 - - - - 
20180916B CHIME 6.03 0.63 506.29 417.98 - - - - 
20180916B LOFAR* 7.68 0.81 3009.12 1313.66 - - - - 
20180916B WSRT.Apertif 7.59 0.78 2664.54 1338.21 - - - - 
20181030A CHIME* 5.15 0.85 247.60 81.49 - - - - 
20190520B FAST 7.70 2.30 681.24 1312.10 -3.22 0.62 0.05 0.05 
20200120E Effelsberg 3.42 1.33 74.35 37.59 0.22 0.63 1.53 9.20 
20201124A CHIME* 5.13 0.63 205.46 185.05 -0.40 0.90 1.00 0.33 
20201124A FAST 4.46 0.96 137.67 111.52 -3.16 0.72 0.06 0.04 
20201124A uGMRT 5.12 0.99 274.93 283.16 - - - - 


NoTE— * The sample size is small. 
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Figure 14. Correlation between various parameters concerning the waiting time for repeating FRBs. u and o are two parameters 


derived by fitting the waiting time through a lognormal distribution for each FRB, characterizing the typical waiting time and 


the dispersion range, respectively. ‘The red dots represent the high-value waiting time while the blue dots are the low-value 
waiting time. Average stands for the mean value. Peak represents the value observed most frequently. The FRBs annotated 


in black were not used for parametric analysis due to the small sample size, but are also shown here for completeness. 


Figure 18 with Figure 15, it could be seen that the amplitude of fluence variation over time is actually very small. 
Meanwhile, the distribution of the fluence is well fitted by a lognormal function. Like the peak flux in Section (3.6), 
for FRB 20121102A, the number of low fluence bursts is slightly higher than the fitting curve for both the Arecibo 
observations (Figure 18: Row 1, Column 2) and the FAST observations (Figure 18: Row 1, Column 3). The large 
number of latest bursts detected from FRB 20201124A by FAST (Zhang et al. 2022b) will help to clarify whether 
other sources have a similar behavior at the low fluence end in this aspect. 


A lognormal function is used to fit the distribution of fluence. The derived parameters are presented in Table 8. Here, 


LLFluence 18 à measure of the typical fluence. It is plotted versus the redshift and SFR for repeating FRBs in Figure 
19. We see that Upriuence 18 negatively correlated with redshift, but positively correlated with SFR. Such correlations 
are likely due to the evolution of source features at high distance scales, which may affect the emission process. Note 
that the data point corresponding to uGMRT detection of FRB 20201124A seems to be an outlier in both panels of 
Figure 19. Its value of HFluence is higher, which may be due to the relatively higher detection threshold of uGMRT 
(Marthi et al. 2022). 
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Figure 15. The peak flux plotted versus the MJD time for repeating FRBs. The inset shows the observational data points of 


different telescopes separately. The horizontal dashed line corresponds to the average of peak flux for each source. The solid line 


is the best linear fit to the data points. The distributions of the peak flux observed by different telescopes are also illustrated 


via histogram for each source. 
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Figure 16. UPeakFlux plotted versus the redshift for repeating FRBs. Here, UPeakFlux is derived by fitting the distribution of 
peak flux with a lognormal function, which is a measure of the typical peak flux for the source. FRBs 20180916B and 20180301A 


are not used for statistical analysis due to the small amount of bursts, but are also shown in the figure for completeness. 


Relations between various statistical parameters concerning the fluence are shown in Figure 20. The dispersion range 


of fluence is negatively correlated with the average fluence, while the average fluence itself is positively correlated with 


the peak. The dispersion range of fluence is negatively correlated with the average fluence (Figure 20: Row 1, Column 
1), but note that FRB 20190520B is still an obvious outlier in the plot. 
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Figure 17. Correlation between various statistical parameters derived by fitting the peak flux distribution with a lognormal 


function for repeating FRBs. u and o characterize the typical peak flux and the dispersion range, respectively. Average stands 


for the mean value of peak flux. Peak represents the value observed most frequently. The FRBs annotated in black are not 


used for statistical analysis due to the small amount of bursts, but are also shown in the figure for completeness. 


Table 7. Parameters derived by fitting the peak flux distribution with 


a lognormal function for repeating FRBs. 


u and ø characterize the 


typical peak flux and the dispersion range, respectively. Average stands 


for the mean value of peak flux. Peak represents the value observed 


most frequently. 


FRB 


20121102A Arecibo* 
20121102A FAST 

20180301A FAST* 

20180916B CHIME* 
20201124A FAST 

20201124A uGMRT 


m c 
(In(mJy)) (In(mJy)) 
3.77 0.95 
2.97 1.34 
3.27 0.75 
0.00003 0.44 

4.56 1.29 
6.72 0.55 


NoTE— * Sample size is small. 


4. DOUBLE-PARAMETER ANALYSIS 


Average 


(mJy) 


68.18 
47.92 
34.70 


220.64 
961.70 


Peak 


(mJy) 


25.13 
10.12 
7.14 
1.51 
56.81 
706.82 


In this section, we analyze the correlation between two parameters among DM, RM, bandwidth, pulse width, waiting 


time, peak flux and fluence, hoping to put forward further constraints on the triggering and emission mechanism of 


repeating FRBs. 
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Figure 18. The fluence plotted versus time for repeating FRBs. The horizontal dashed line corresponds to the average fluence 
of each FRB. The solid line is the best linear fit to the data points. The inset shows the observational data points from different 


telescopes separately. For each source, the distribution of the fluence is also illustrated via histogram according to the detection 


telescope. 


4.1. Pulse Width, Peak Flux and Fluence 


The pulse width, peak flux and fluence are three basic parameters of FRBs, which are strongly connected with the 
energetics of the burst. Fluence is often calculated by multiplying the width with the peak flux. In Figure 25 FRBs 
are plotted on the peak flux - pulse width plane. We see that is no obvious correlation between these two parameters, 
which means that the peak flux and the width are two independent parameters. An interesting feature is that, for 
each FRB source, the peak flux typically varies in a wide range (spanning for more than two magnitudes), but the 
width only varies in a much narrower range (spanning for less than one magnitude). This is especially clear for the 
two most active sources of FRBs 20121102A and 202011244. In Figure 26, FRBs are plotted on the width - fluence 
plane. There is a weak positive correlation between the width and fluence for each source. It is mainly due to the fact 
that the fluence is calculated by multiplying the width with the peak flux. However, note that Xiao & Dai (2022a) 


suggested an alternative explanation for the weak correlation between pulse width and fluence for FRB 201211024, 
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Figure 19. Correlation between redshift, SFR and the typical fluence (uriuence) for repeating FRB sources. Panel (a): Wriuence 


versus redshift. Panel (b): uriuence versus SFR. Uriuence is derived by fitting the distribution of fluence with a lognormal function 


for each FRB, which represents the typical fluence. FRB 201803014 is not used for statistical analysis due to the small sample 


size, but is also shown in the figure for completeness. 
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Figure 20. Correlation between various statistical parameters derived by fitting the fluence distribution with a lognormal 


function for repeating FRBs. u and o characterize the typical fluence and the dispersion range, respectively. Average stands 


for the mean value of fluence. Peak represents the value observed most frequently. The FRBs annotated in black are not used 


for statistical analysis due to the small amount of bursts, but are also shown in the figure for completeness. 
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Table 8. Parameters derived by fitting the fluence distribution with a lognor- 


mal function for repeating FRBs. p and o characterize the typical peak flux 


and the dispersion range, respectively. Average stands for the mean value of 


peak flux. Peak represents the value observed most frequently. 


FRB m c Average Peak 
(In(Jyms)) (In(Jyms)) (Jyms) (Jyms) 
20121102A Arecibo -2.03 1.42 0.36 0.19 
20121102A FAST -2.51 1.32 0.19 0.05 
20121102A Lovell 0.42 0.93 2.35 0.92 
20180301A FAST* -2.18 0.53 0.13 0.18 
20180916B CHIME* 1.16 1.03 5.41 2.36 
20180916B LOFAR* 4.96 0.51 162.89 109.93 
20180916B WSRT.Apertif 1.94 0.88 10.21 4.80 
20190520B FAST -2.26 0.50 0.12 0.11 
20200120E Effelsberg -1.94 0.60 0.17 0.10 
20201124A FAST -0.25 1.25 1.70 0.34 
20201124A uGMRT 2.69 0.80 20.21 12.28 


NoTE— * Sample size is small. 


namely the bursts may be divided into two subtypes according to their brightness temperature and may be produced 


by different physical processes. FRBs are plotted on the peak flux - fluence plane, in Figure 27. We could see a clear 


positive correlation between the peak flux and the fluence. Again, the reason is that the fluence is a product of the 


peak flux and the width, with the width varying only in a narrow range. 


Correlation between the average pulse width (uwiath), peak flux (UPeakFlux) and fluence (WFiyence) is plotted in 


Figure 21 (Row 1) for repeating FRB sources. These parameters are also plotted versus DM Host and | RM|. There 


is a positive correlation between the average peak flux and pulse width (Figure 21: Row 1, Column 1). However, this 


correlation may be due to the selection effect of observations, since a wider pulse and a higher peak flux will make the 


FRB more easy to be detected. UPeakFluzr aNd UMriuence are tightly correlated as shown in Figure 21 (Row 1, Column 
2), while the correlation between pwidth and HFluence i8 very weak (Figure 21: Row 1, Column 3). Again it is mainly 
due to the fact that UpeakFlux Varies in a wide range, but wiadtn 1s distributed in a much narrower range. As a result, 
the influence of peak flux on fluence is more significant than that of pulse width. Additionally, UWiatn is affected 
by cosmological expansion and propagation effects (see Section 3.4), which further interferes the correlation between 


UWidth aNd UFiuence. Note that the linear correlation between Lwigth, MPeakFlux aNd HFluence actually corresponds to 


a power-law function among the average of width, peak flux and fluence, which is roughly consistent with the results 
of Xiao & Dai (2022b). 

Figure 21 (Row 2, Column 1) shows that uwiatn is positively correlated with DM Host for most FRB sources. This 
is not strange, since a larger DM tends to broaden the pulse width. However, note that FRB 20190520B markedly 


deviates from this general trend, the reason of which is still unknown. 


4.2. Bandwidth, Pulse Width and Waiting Time 


The bandwidth, pulse width and waiting time are all influenced by cosmological expansion, and might also be related 


to potential redshift-dependent evolution of the central engine (see Section 3.3, 3.4, 3.5). Correlations between the 
average bandwidth (WBandwidtn), pulse width (uwiatn) and waiting time (UwaitingTime(High—Value)) are illustrated in 
Figure 22 for repeating FRB sources. These parameters are also plotted versus DM Host and | RM |. From Row 1 of this 
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Figure 21. Correlation between the average pulse width (wwiatn), peak flux (UPeakFiux) and fluence (UFiuence) for repeating 
FRB sources. These parameters are also plotted versus DM gos; and |RM|. Here, DM Host represents the DM component 
contributed by the host galaxy and the local circum-burst environment. The FRBs annotated in black were not used for the 


double-parameter analysis due to the small sample size, but are also shown in the figure for completeness. 


figure, we see that Bandwidth; HWidth. aNd LwaitingTime(High—Value) are positively correlated with each other. Note 
that the Bandwidth - HWiath correlation (Figure 22: Row 1, Column 1) and the Bandwidth - HWaitingTime(High-Value) 
correlation (Figure 22: Row 1, Column 2) are both very tight, except that the Arecibo data point of FRB 20121102A 
seems to be an obvious outlier in the two panels. The deviation of FRB 20121102A (Arecibo observations) may be due 
to the relatively wider work bandwidth of Arecibo, which is 1.15 — 1.73 GHz). As a comparison, the work bandwidth 
of Effelsberg is 1.2 — 1.6 GHz and the band of FAST is 1.05 — 1.45 GHz (see Section 3.3). However, the uwiatn - 
UW aitingTime(High—Value) Correlation is very dispersive as compared with the above two correlations. It is not clear 
whether such a dispersion is due to the increased sample size or not. In the future, more observational data points are 
needed to further clarify this issue. In Figure 22 (Row 2, Column 1), there is a weak correlation between Bandwidth 


and DM Host. In other panels of Figure 22, no clear correlation is observed. 


4.3. Bandwidth, Waiting Time and Peak Flux, Fluence 


Correlations between the average bandwidth (upanawiatn), waiting time (Lwaitingrime(High-Value)); peak flux 


(LPeakFluz), and fluence (UFluence) are illustrated in Figure 23 for repeating FRB sources. In Panels (a) and (b), 
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Figure 22. Correlation between the average bandwidth (UBandwidth), pulse width (uwiath) and waiting time 
(LW aitingTime(High—Value)) for repeating FRB sources. These parameters are also plotted versus DM Host and |RM|. Here, 
DM Host represents the DM component contributed by the host galaxy and the local circum-burst environment. The FRBs 


annotated in black were not used for the double-parameter analysis due to the small sample size, but are also shown in the 


figure for completeness. 


since the observational data points are too few, no firm conclusion could be drawn on the correlation between the 
corresponding parameter pairs. In Panel (c) there seems to be a positive correlation between Mw aitingTime(High—Value) 
and peak flux (HPeakFlux), especially when the data point corresponding to FAST observations of FRB 201803014 is 


expelled due to its small sample size. However, be aware that the available data points are still very limited in this 


panel. A weak positive correlation could also be interestingly observed between Mw aitingTime(High—Value) ANd HFluence 
in Panel (d). It means that for different repeaters, when the high-value waiting time is larger, the energetics of bursts 


is also larger, which may be connected with the triggering mechanism of repeating FRBs. 


4.4. Dispersion Measure and Rotation Measure 


DM and RM can reflect the physical conditions of the surrounding media around the source region. For repeating 
FRBs with a measured redshift, we can estimate DM Host based on the Macquart DM rGm — z relation and the MW 
free electron models such as NE2001 or YMW16 (see Section 3.1). Note that DM Host refers to the DM contributed 
by the host galaxy and the local circum-burst environment of the FRB source. 

Correlation between |RM| and DM Host are plotted in Figure 24 for repeating FRBs. It is striking that these two 
parameters are tightly correlated and can be well fitted by a linear function. Such a linear correlation strongly indicates 
that the magnetic field in the local environment does not vary too much for different FRB sources. However, note 
that FRB 20121102A deviates from the best fit line markedly. FAST recently detected a series of bursts from FRB 
20121102A again, with the DM measured to be significantly smaller (Wang et al. 2022b), which makes the deviation 
even more significant. It seems that this repeater may have a much stronger magnetic field. 

The linear correlation between IRM | and DM Host can potentially help us to estimate the distance of those FRBs 
whose redshift is not measured. For this purpose, we first need to measure the | RM |. Then the |RM| - DM gos; can be 
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Figure 23. Correlation between the average bandwidth (upanawiatn), waiting time (IwaitingTime(High—Value)), peak flux 
(UPeakFlux), and fluence (UFiuence) for repeating FRB sources. The FRBs annotated in black were not used for the double- 


parameter analysis due to the small sample size, but are also shown in the figure for completeness. 


applied to estimate the DM Host. Finally, we can calculate DM rGm and use it to estimate the distance. However, this 

approach is mainly limited to repeating FRBs, because measuring | RM | requires observations over long-time scales. 
In our study, we have also explored potential correlations between many other parameter pairs, but fail to find any 

clear correlation among them. More double-parameter distributions of FRBs where no correlation is observed are 


illustrated in the Appendix Figures of 28 — 51. 


5. CONCLUSIONS AND DISCUSSION 


In this study, a sample of 16 FRB repeaters is established, with at least five bursts being detected from each 
source. Among them, four repeaters, i.e. FRBs 20121102A, 20180916B, 20190520B and 20201124A, are the most 
active sources. More than 5,000 bursts are included, which are observed by 11 different telescopes. Data of many 
key parameters are collected from the literature, including the trigger time, pulse width, dispersion measure (DM), 
Faraday rotation measure (RM), bandwidth, waiting time, peak flux and fluence. A thorough statistical investigation 
is performed based on the available data and the behaviors of different repeaters are compared. The DM of FRB 
201211024 is found to increase continuously on a long timescale. For most repeating sources, DM typically varies in a 


very narrow range of +3 cm ^? pc. However, FRB 20190520B is found to have a large variation of +12 cm ^? pc. The 
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Figure 24. Correlation between |RM| and DM mos: for repeating FRBs. DM Host represents the DM component contributed 


by the host galaxy and the local circum-burst environment. 


Macquart relation (Macquart et al. 2020b), i.e. the tight positive correlation between DM ram and z, is confirmed 
with this expanded sample. Furthermore, DM is also found to be linearly correlated with SFR, but the absolute value 


of RM correlates negatively with SFR. RM evolves with time in a much more chaos behavior in different repeaters. 
The pulse width generally does not evolve with time. The waiting time shows a bimodal distribution in most repeaters. 
The absolute RM is also found to be linearly dependent on DMhost (but note that FRB 20121102A is an obvious 
outlier), which possibly hints the existence of a universal magnetic field in different sources. 

Contrary to the increasing tendency of DM evolution in FRB 20121102A, the DM of FRB 20190520B seems to be 
decreasing over a long time, although the decreasing tendency is rather weak. Such a decreasing, together with the 
relatively wide variation range of DM mentioned just above, indicates that it is quite special as compared with other 
repeaters. In fact, FRB 20190520B deviates from other sources in many plots in Figure 2. It also has an extremely 
large DMhost, which is ~ 1208 cm ? pc. It even deviates from the Macquart relation seriously (see Figure 2). The 
intrinsic reasons that lead to these special features are still unknown. Note that the number of bursts detected from 
FRB 20190520B is still not too large, i.e. about 150 in total. More observations on this repeater in the future will be 
helpful for understanding its nature. 

Only two repeaters have well measured RM values in a long period, i.e. FRBs 20201124A and 20180916B. The RM 
of FRB 20201124A evolves in a very special pattern over a time range of ~ 54 days (Figure 2), with an amplitude of 
~ 540 rad m~?. On the contrary, the RM evolution is quite simple for FRB 20180916B, which only monotonously 
increases by a small amplitude of ~ 50 rad m^? over a long period of ~ 1000 days. It clearly indicates that there 


might be more than one kind of triggering mechanisms even for repeating FRBs. 


For the repeaters in our sample, the average bandwidth is found to be closely related to the average pulse width 


and the average waiting time. However, no correlation exists between the average pulse width and the average waiting 


time, which means that the cosmological time dilation effect is not significant for the sample. The correlation between 


the average bandwidth and the average pulse width/average waiting time is then likely due to the evolution of the 


source characteristics at different redshifts. 
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APPENDIX 


In addition to the above analyses, we have also explored a few other double-parameter distributions of repeating 


FRB sources. But generally, no clear correlation is observed. Here, we present these additional double-parameter plots 


as an Appendix. 
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Figure 25. The distribution of FRBs on the peak flux - pulse width plane. Bursts from each source are plotted according 
to the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 26. The distribution of FRBs on the width - fluence plane. Bursts from each source are plotted according to the detection 


telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, respectively. 
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Figure 27. The distribution of FRBs on the peak flux - fluence plane. Bursts from each source are plotted according to 
the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 


38 HU ET AL. 


F T lp =] E 
4 FRB20121102A 102 i : FRB20121102A m "j FRB20121102A x FRB20180301A 
Z Arecibo g FAST £ F 1 Lovell £ FAST 
S 2 S io (1.05-1.45GHz) S (1.312-1.712GHz) 5 (1.05-1.45GHz) 
o . o === Average o ==» Average o ==» Average 
o ur: o u o Bursts o Bursts 
^ 10? Eb : = 10° = T——T 4 0 
566 I 1 570 } ° I 525.0 I 
I , 570+ I 4 I 
201 E DC ME Pm mar ! 522.5 - l 
E l — 565F e — 568 I + TERE I 
e 562 LI g 7 g I 9. I 
E: I «- 566r I iae Dried, lan Pa x PME EID 
! 560 I D 1 I D [| 
E I E seo k E seat I 5 517.5 I 
5^ A au uo D aer a r€— E = I = e | = I 
= 1 = 1 A ser sT Š 515.0 : 
556 - l 555 F ! l I 
[| [| 560+ [| [| 
512.5 
554 I I I I 
[| [| 558 [| [| 
i 1 i 550k. i sul sul L 1 i 1 fer 510.0 e 1 
10! 0 2 4 10° 10! 10 10! 1o 109 10? 10! 10! 0 5 
Width (ms) Counts Width (ms) Counts Width (ms) Counts Width (ms) Counts 
T T p T T T 10? T T T 
FRB20180916B FRB20190520B FRB20201124A 
2 24 CHIME 2 10! FAST 2 FAST 
S (0.4-0.8GHz) E d (1.05-1.45GHz) S 10 (1.05-1.45GHz) 
o == Average o == » Average o == Average 
o Bursts o Bursts o Bursts 
OF 4 - + T 109 E pu. + + — - + + 10° T mr y 
I x I 
349.8 - l 1 422.5 
I I 
^ 349.6 F l l — 420.0 
o j [| [| o 
o I I a 
m m 
|. 349.4 l si P 4175 
E H----------- 4 E 
= l E 415.0 
= 349.2 F l = : 
a i a 
412.5 
349.0 - l 
l 
I 410.0 
348.8f e , L 
10° 10! 1 10! 10? 
Width (ms) Counts Width (ms) Counts Width (ms) Counts 


Figure 28. The distribution of FRBs on the DM - pulse width plane. Bursts from each source are plotted according to 
the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 29. The distribution of observed bursts on the RM - pulse width plane, for FRB 20180301A and FRB 20201124A. 
All the bursts are detected by FAST. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 30. The distribution of FRBs on the DM - peak flux plane. Bursts from each source are plotted according to 


the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 31. The distribution of observed bursts on the RM - peak flux plane, for FRB 20180301A and FRB 20201124A. All 


the bursts are detected by FAST. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 32. The distribution of FRBs on the DM - fluence plane. Bursts from each source are plotted according to the detection 


telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, respectively. 
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Figure 33. The distribution of observed bursts on the RM - fluence plane for FRB 20201124A. All the bursts are detected by 


FAST. The horizontal and vertical dashed lines mark the average values of the two parameters, respectively. 
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Figure 34. The distribution of FRBs on the bandwith - pulse width plane. Bursts from each source are plotted according 
to the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 35. The distribution of FRBs on the bandwith - waiting time plane. Bursts from each source are plotted according 
to the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 36. The distribution of FRBs on the pulse width - waiting time plane. Bursts from each source are plotted according 


to the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 37. The distribution of observed bursts on the DM - bandwidth plane, for FRB 20121102A and FRB 20201124A. 


All the bursts are detected by FAST. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 38. The distribution of observed bursts on the RM - bandwith plane for FRB 20201124A 


by FAST. The horizontal and vertical dashed lines mark the average values of the two parameters, 


. All the bursts are detected 


respectively. 
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Figure 39. The distribution of FRBs on the DM - waiting time plane. Bursts from each source are plotted according to 


the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 40. The distribution of observed bursts on the RM - waiting time plane, for FRB 20180301A and FRB 20201124A. 


All the bursts are detected by FAST. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 41. The distribution of observed bursts on the bandwidth bound - pulse width plane for FRB 20121102A. All the 
bursts are detected by Arecibo. fHigh and frow stand for the upper and lower bounds of observed bandwidth, respectively. The 


horizontal and vertical dashed lines mark the average values of the two parameters. 
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Figure 42. The distribution of observed bursts on the bandwidth bound - waiting time plane for FRB 20121102A. All the 
bursts are detected by Arecibo. fHigh and frow stand for the upper and lower bounds of observed bandwidth, respectively. The 


horizontal and vertical dashed lines mark the average values of the two parameters. 
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Figure 43. The observed center frequency potted versus pulse width for FRB 20121102A 
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. All the bursts are detected by 


Arecibo. The horizontal and vertical dashed lines mark the average values of the two parameters, respectively. 
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Figure 44. The distribution of FRBs on the bandwidth - peak flux plane. Bursts from each source are plotted according to 


the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 45. The distribution of FRBs on the bandwidth - fluence plane. Bursts from each source are plotted according to 


the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 46. The distribution of FRBs on the peak flux - waiting time plane. Bursts from each source are plotted according 


to the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 


46 HU ET AL. 


FRB20121102A FRB20180916B 


FRB20121102A 10? FRB20121102A 
FAST 


g Arecibo 2 2 Lovell 8 WSRT.Apertif 
5 (1.15-1.73GHz) 5 (1.05-1.45GHz) 5 (1.312-1.712GHz) E ] (1.22-1.52GHz) 
o o 1 Bursts o o 1 Bursts 
O O 10 = « Average o = » Average o = + Average 
10° 
10? I 
l 
I 
= = ax =~ I 
u u I] Uu ;| 
E E E E t 
3 3 3 "Ep a ee quM S a ne oS 
[^] [^] (MUI cce on ie mda ie ne ba jer ee dl aite. [^] [| 
o o o o 
c c c c k 
[7] [7] [7] [7] I 
= =] = E 
iE u i 10° u 1 
i 
I 
l 
10° 
107? 107! 109 10! 10? 10 10° 10! 10? 10-3 107! 10! 10 10 10! 109 10? 103 10° 10? 103 10* 10° 10! 
Waiting Time (s) Counts Waiting Time (s) Counts Waiting Time (s) Counts Waiting Time (s) Counts 
pee i i j i "] FRB201905208 10: m f zn FRB20200120E ioi FRB20201124A 102 FRB20201124A 
2 aol J FAST 2 Effelsberg P uGMRT E 
5 E J (1.05-1.45GHz) S (1.2-1.6GHz) S (0.55-0.75GHz) 10 (1.05-1.45GHz) 
o Bursts o ursts o Bursts Bursts 
eo =» Average o * Average o = » Average 
+ + — 109 —Á: 10° Ẹ -+ TU 109 
I 4 10? [| 
| l 
| I 
R = l CMM I 10! 
N u l u I 
E E i E i 
>» ERe———-———— ooo omms mc m m DD n m > > 
> 10-1 i^ D — þe —————— LG „Pab am m a a a a => qBR—————————— aa ÉBP——-—T-—— | | f[pr-— 
v v l v l 109 
o o o 
5 T 10-1 © 10 - 
5 as r Du l 
ma [a p us f 
I i I 1071 
I I 
I | I 
1 1 n Luul È 1 fi 1 1 L 1 1 1 | en suit sul L "TERI I 1 eo soit 1 mI 1 " f ut 1 sul mm 
107? 107! 10° 10! 10? 103 10! 10° 10! 10? 103 10! 10° 10! 10? 103 10! 10-2 107! 109 10! 10? 10? 109 10! 10? 
Waiting Time (s) Counts Waiting Time (s) Counts Waiting Time (s) Counts Waiting Time (s) Counts 


Figure 47. The distribution of FRBs on the fluence - waiting time plane. Bursts from each source are plotted according to 
the detection telescope separately. The horizontal and vertical dashed lines mark the average values of the two parameters, 


respectively. 
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Figure 48. The distribution of observed bursts on the bandwidth bound - fluence plane for FRB 20121102A. All the bursts are 
detected by Arecibo. fHigh and frow stand for the upper and lower bounds of observed bandwidth, respectively. The horizontal 


and vertical dashed lines mark the average values of the two parameters. 
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Figure 49. The observed center frequency potted versus flux density for FRB 20121102A. All the bursts are detected by 


Arecibo. The horizontal and vertical dashed lines mark the average values of the two parameters, respectively. 
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Figure 50. The observed center frequency potted versus fluence for FRB 20121102A. All the bursts are detected by Arecibo. 


The horizontal and vertical dashed lines mark the average values of the two parameters, respectively. 
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Figure 51. The distribution of observed bursts on the DM - RM plane for FRB 20180301A and FRB 20201124A. All the bursts 


are detected by FAST. The horizontal and vertical dashed lines mark the average values of the two parameters, respectively. 


